Astronomy & Astrophysics manuscript no. specevol 


February 5, 2008 


(DOI: will be inserted by hand later) 





Long Term Variability of Cyg X-1 

IV. Spectral Evolution 1999-2004 

J. Wilms\ M.A. Nowak^, K. Pottschmidt^ G.G. Pooley"^, and S. Fritz^ 

' Department of Physics, University of Warwick, Coventry, CV7 lAL, UK 

2 MIT-CXC, NE80-6077, 77 Massachusetts Ave., Cambridge, MA 02139, USA 

^ Center for Astrophysics and Space Sciences, University of California at San Diego, La Jolla, CA 92093-0424, USA 

'* MuUard Radio Astronomy Observatory, Cavendish Laboratory, Madingley Road, Cambridge CBS OHE, UK 

^ Institut fiir Astronomic und Astrophysik, Universitat Tubingen, Sand 1, 72076 Tubingen, Germany 

Received <date> / Accepted <date> 

Abstract. Continuing the observational campaign initiated by our group, we present the long term spectral evolution of the 
Galactic black hole candidate Cygnus X-1 in the X-rays and at 15 GHz. We present ~200 pointed observations taken between 
early 1999 and late 2004 with the Rossi X-ray Timing Explorer and the Ryle radio telescope. The X-ray spectra are remarkably 
well descri bed by a simple b roken power law spectrum with an exponential cutoff. Physically motivated Comptonization mod- 
els, e.g., by'Titarchuk' ^1994', compTT) and bv lCoppil il999^. eqpair), can reproduce this simplicity; however, the success of the 
phenomenological broken power law models cautions against "overparameterizing" the more physical models. Broken power 
law models reveal a significant linear correlation between the photon index of the lower energy power law and the hardening 
of the power law at ~10keV. This phenomenological soft/hard power law correlation is partly attributable to correlations of 
broad band continuum components, rather than being dominated by the weak hardness/reflection fraction correlation present 
in the Comptonization model. Specifically, the Comptonization models show that the bolometric flux of a soft excess (e.g., 
disk component) is strongly correlated with the compactness ratio of the Comptonizing medium, with Ldisk (4/4)"""". Over 
the course of our campaign, Cyg X- 1 transited several times into the soft state, and exhibited a large number of "failed state 
transitions". The fraction of the time spent in such low radio emission/soft X-ray spectral states has increased from ~10% in 
1996-2000 to ~34% since early 2000. We find that radio flares typically occur during state transitions and failed state transi- 
tions (at /"h/^s ~ 3), and that there is a strong correlation between the 10-50 keV X-ray flux and the radio luminosity of the 
source. We demonstrate that rather than there being distinctly separated states, in contrast to the timing properties the spectrum 
of Cyg X-1 shows variations between extremes of properties, with clear cut examples of spectra at every intermediate point in 
the observed spectral correlations. 
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1. Introduction 

One of the major results of X-ray and y-ray astronomy since the 
discovery of Galactic black hole candidates (BHCs) more than 
thirty years ago has been the realization of the complexity of 
the spectral variability of these sources on time scales from mil- 
liseconds to decades. It is generally believed that this variabil- 
ity can provide clues as to the basic physical processes in X-ray 
binaries. Importantly, the broad-band spectra and timing behav- 
ior of accreting stellar mass black holes in X-ray binaries and 
of supermassive black holes in the centers of Ac t ive Galactic 
Nucle i (AGN) are very similar dUttlev et a l.' 2002-. Uttlev etaP 
l2005b . Due to the scaling of the dynamical time scales in these 
systems with the mass of the central object, studying Galactic 
and extragalactic systems is complementary. Galactic sources 
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allow us to see how the accretion process evolves on very long 
time scales, while AGN provide "snapshot" observations of the 
accretion process, over characteristic dynamical time scales, 
with a signal to noise ratio that is not reachable for Galactic 
sources. 

Observations over the last three decades have shown 
that Galactic BHC exhibit distinct and very characteristic 
states with different spectral shape and varia bility behavior, 

with AGN seemingly following similar trends ( IMcHardv et alJ 

I I I t I ° ° 
2003 Jester 2005, and therein). The sta te of a BHC depends on 

not yet fully understood parameters tHoman et al...2001.) . al- 
though it is believed that the mass accretion rate, M, and there- 
fore tiieJuminosit^_rf_fli^^ system, plays a vital role. 
See lMcClintock & Remiii^ 02OO3.) and .BellonL (.2004.) for an 
extensive discussion of these issues. 
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For low accretion luminosities in the "hard state" of 
Galactic black holes, the X-ray spectrum can be described by 
a hard power law with photon index T ~ 1.7 and an expo- 
nential cutoff at ~150keV. At low energies, most BHC show 
some kind of soft excess with a characteristic temperature of 
a few lOOeV. This baseline continuum is modified by a Fe 
Ka emission line at ~6.4keV and by reflection features, in- 
dicating the close proximity of the source of hard X-rays and 
relatively cold material. The hard state is furthermore charac- 
terized by strong X-ray variability of >10% rms and by the 
presence of radio emission. In systems where the radio sources 
have been resolved, this radio emission has been shown to 
originate in an outflow that is consistent with being mildly 
relativistic (v ~ 0.6c) in the hard state (iStirling et all 1200 it 
iGallo et alJ2003l) . For Cyg X-1, this outflow has recently been 
shown to have a time averaged kinetic power of ~10^* erg s ' 
to 10^^ erg s~^ i.e., a significa nt fraction of the system's X-ray 
luminosity l OaUo et alJ2005h . 

In contrast to the "hard state", the (typically) higher lu- 
minosity "soft state" exhibits a soft X-ray spectrum that can 
be well described by thermal emission from a standard accre- 
tion disk. Where a hard spectral component is detected, it does 
not sho w any appreciable curvature, sometimes up to the MeV 
regime (IPro ve et al.lfl998l: iMcConnell et alJ^OOoi l2002i) . The 
X-ray variability during the s oft state is weak (< 6% rms) and 
no radio emission is detected jFender et alJ 19991) . 

While there is general agreement that the soft excess is due 
to emission from an accretion disk, the interpretation of the 
hard spectral component and its relationship to the reprocess- 
ing features is still debated. The canonical interpretation of the 
hard state X-ray spectrum, first proposed in the 1970s (e.g.. 



hard state A-ray spectrum, iirst proposed m the 19/Us (e.g., 
i Thorne & Pricell975lShapiro et alJl976HSunvaev & Triimpei 



|1222L ™d Isiter elaborated upon, e.g., bv iHaardt & Marasch: 
lll99ll) and Haardt & Maraschi ( 1993), is that the hard spectral 
component is caused by Comptonization, where soft X-rays 
from the inner disk are Compton upscattered by hot electrons 
(^re ~ lOOkeV) in a predominantly thermal electron gas, of- 
ten called the "accretion disk corona" (ADC). The spectrum 
emerging from the corona has the proper power law plus expo- 
nential cutoff shape. To explain the pr esence of the ADC, mag- 
netohydrodynamical instabilities (e.g.. Balbus & Hawlev^ l998l 
and references therein), which only work efficiently at lower M, 
have been invoked. 

Observations limit the covering factor of the ADC with 
respect to the source of soft X-rays to «:1. Otherwise 
the ADC would fail to reach the high coronal tempera- 
tures, inferred from the energy of the exponential cutoff, 
due to efficient Co mpton cooHng (e.g.. iHaardt et alJ Il997t 
iDove et al]ll997alb[ and therein). Geometric models invoked 
to achiev e smaller ADC covering factors in clude patchy 
coronae (*Stern e t ^ all fl99Sl magnetic flares d Beloborodovl 
Il999t [Poutanen & Fabianlll9 99). advection dominated accre- 
tion flows (Esin et al. 1998) , ADCs with non-sta t ic or outflow- 
ing c oronae (.Beloborodov & Illarionovl l200lt iMalzac et al.l 
12001 1) , and accretion flows with an inner geometrically thick 
ADC and an outer geometrically thin and optically thick accre- 
tion disk (.Dove et al.. 1997 a b). In all of these models, a frac- 
tion of the hard X-rays is scattered back towards the disk, giv- 



ing rise to fluo rescence Fe Ka emission and the Compton re- 
flectio n hump jLi ^htman & Rvbic'killl979t iLi ghtman & Whitel 
'1988'). Comparisons between observations and theoretical 
model spectra for all of these different geometries have been 
successful. 

With the realization of the importance of jets and their as- 
sociation with the hard state, the "Comptonization paradigm" 
has recently been challenged. The discovery of a tight correla- 
tion between the radio and the X-ray emission on time scales of 
days (Hannikainen et al. 1998; Poolev et al. 1999; Corbel et alj 
I2OOOI) suggests a clos e coupling between t he X-ray and the 
radio emitting media jMarkoff et al.l 12003^ iHeinz & Sunvaev| 
2003). In addition, in a direct comparison between data and 
theory, jet models capable of reproducing the radio-X-ray cor- 
relations have been recently shown to describe the X-ray spec- 
trum with a precisio n comparable to Comptonization models 
llMarkofl^etalJl200.5l) . The origin of the X-rays in these mod- 
els is due to a combination of synchrotron emission from the 
jet and synchrotron self Compton emission from the jet base. 
We note, however, that a small separation of the emission re- 
gions cannot be ruled out from ioint radio-X-ray timing tim- 
ing arguments (Gleissneret al. 2004a), furthermore, accretion 
models postulating a strong coupling between a Comptonizing 
medium and the jet are also able to explain the radio-X-ray cor- 
relation without postulating a jet origin for the X-rays (Meier] 
l200lUMerloni & Fabianl2002l) . 

In conclusion, all existing models suggest a complex inter- 
play between the energetics and emission from the accretion 
disk, Comptonizing plasma, and radio jet. This interplay can 
be disturbed by changes in some external parameter, such as 
M. These changes may lead to state changes as well as to subtle 
changes in the overall so urce properties, such as the so called 
"failed state transitions" jPottschmidt et aljl200(i) . An empiri- 
cal understanding of the overall accretion process is therefore 
difficult to gain from only a few single and isolated observa- 
tions, although these are very important for determining broad 
band spectr a or for measuring spectra with high energy reso- 
lution (e.g.. lMilleretalJl2002 ). Monitoring campaigns cover- 
ing the characteristic time scales of these spectral changes are 
required. We note that the necessity of monitoring campaigns 
follows also from hysteresis effects seen in many sources, 
i.e., the earlier source history is important for determining the 
source properties in a given observation ( Miyamoto et al. 19951 
iNowak et a l. 2002: Maccarone & C oppi 2003). 

Monitoring campaigns using large effective area and pro- 
portional counter energy resolution became possible in the 
1990s with the Rossi X-ray Timing Explorer (RXTE). Since 
then, monitoring campaigns using RXTE's pointed ins tru- 
ments , the Proportional Counter Array, (PCA; Jahoda e t alJ 
I1997') an d the High Energy X -ray Timing Experiment, 
(HEXTE; iRothschild etsJl 1 19981) . have been performed for 
virtually all known persistent and transient Galactic black 
hole candidates. Among many others, examples for s uch cam- 
paign s are those on LMC X-1 and LM C X-3 (JVowak et alJ 
200 it IWilms et al.| 1200 ll). GX 339-4 ijWilms et alJ Il999l: 
Nowaketal. 1999, 2002; Kon g et alJ 12002'; Homan et al.j 
2005: Belloni et al.. .20051 V1408 Aal (.Nowak & Wilmsi 
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jKalemci et alJl2003l) . 4 U 1630-47 jTomsick & Kaarel2000) 
and XTE J 1550-564 JKalemci et alJ 1200 ll: LSobczak et alJ 
I2OOOI) . 

In this paper, we present results from one of the longest of 
these monitoring campaigns, which was initiated by the authors 
in 1996 to be performed on the prototypical hard state black 
hole candidate (BHC) Cygnus X-1. Our campaign consists of 
biweekly, ~5 ksec long RXTE pointings and simultaneous ra- 
dio observations at 15 GHz (A -2 cm) with the Mullard Radio 
Astronomy Observatory's Ryle telescope in Cambridge, UK. 
Earlier papers in this series concentrate d on the evolution of X- 
ray time lags ('Pottschmidt et alJl2000i) . of the power spectrum 
(Pottschmidt et al. 2003b, hereafter paper i), the linear relation- 
ship bet ween the short term ro ot mean square variability and 
the flux (i Oleissner et alJ2004 bl. hereafter paper 11), and correla- 
tions between the soft X-rays and the radio flux ( Gleissner et al. 
l2004al hereafter paper m). Amongst others, archival data from 
the campaign have also been used in the discovery of giant 
flares from Cyg X-1 (Gierliiiski & Zdziarski 2003), a subset 
of the archival observations was also u sed in studies of Cyg 
X-l's spectral and temporal variability (" Axelsson et"al]l2005[ 
[ibragimov et al. 2005). Here, we discuss the results obtained on 
the spectral evolution of Cyg X- 1 throughout the RXTE mon- 
itoring campaign, starting in Sect. |2] with a description of the 
changes of our data reduction procedure with respect to pa- 
pers I through III. In Sect.|3lwe give the results of modeling the 
broad band 3-200 keV X-ray spectra using empirical spectral 
models as well as advanced Comptonization models and con- 
sider the evidence for interaction between the X-ray and the 
radio emission. We summarize our results in Sect. |3 We will 
discuss the joint spectral-timing behavior of Cyg X-1 based on 
our analyses in a future paper. 

2. Data Reduction 

2.1. Overview 

We have already given an extensive overview of the RXTE ob- 
serving strategy and our data extraction in papers i through m. 
Here, we only highlight the most important points pertaining 
to the X-ray spectral analysis, mainly centered on the RXTE 
calibration. 

In this paper we consider 202 observations with a typical 
good time of 3-7 ksec, for a total good time of ~989ksec. 
We use data from both the PCA and the HEXTE. Spectra 
are extra cted with HEA SOFT 5.3.1 and then fit with XSPEC 
11.3.1w (Arnaud 1996). For the spectral analysis, we gener- 
ated 3-25 keV PCA spectra from the top Xenon layer data and 
analyze spectra taken with different numbers of proportional 
counter units (PCUs) separately. This choice was made to al- 
low the comparability with the results of the timing analysis. 
For the HEXTE, data from 18-120 keV were considered. 

As shown in Appendix IaI compared to earlier versions of 
the response matrices, the PCA is now in much better agree- 
ment with the HEXTE and also with other missions, in both, 
flux and spectral slope, with the remaining uncertainty being 
taken into account by a multiplicative constant that is nor- 
mally very close to unity. Our reevaluation of the PCA cali- 



bration also shows that remaining calibration uncertainties can 
be modeled to first order by adding a systematic uncertainty 
of 0.5% to the data, although even with this choice some sig- 
nificant calibration effects are still visible in the residuals (see 
AppendixlAt. As shown in AppendixiBl however, this system- 
atic error leads to a dramatic overestimation of the range of the 
confidence intervals for the fit parameters. Thus, although we 
compute confidence intervals at the 90% level for one interest- 
ing parameter, we do not show error bars in the figures. 

The good time of the observations is mainly defined by ex- 
cluding times of high PCA background. Due to the require- 
ments of X-ray timing analysis, early in the campaign very 
conservative constraints were chosen, i.e., data taken within 
30 minutes of passages through the South Atlantic Anomaly 
(SAA) and during times of high particle background were dis- 
carded. While in principle improved PCA background models 
now allow for good results even closer to the SAA passages, 
for consistency reasons and in order to facilitate the compari- 
son of spectral parameters to the results of papers i through iii, 
the previous limits were retained. 

Electronic tables in the Flexible Image Transport System 
(FITS) format containing the results of the spectral fitting and 
the confidence intervals for all best fit parameters are available 
accompanying the electronic version of this paper'. In our dis- 
cussion of peculiar observations below, we identify observa- 
tions in a way that makes it easy to find them in these FITS 
files and in the RXTE archive. The abbreviated syntax used is 
of the type Pxxxxx/XX where xxxxx is the RXTE proposal ID 
and where XX is the number of the observation within each 
proposal. We also give the date and time of the observation to 
the closest hour and, in cases where an observation resulted in 
more than one spectrum because different detector combina- 
tions were used, we identify those detectors that were off. 

2.2. Fitting Strategy 

For larger data sets, the CPU time needed to perform spectral 
fitting can be considerable and therefore the choice of start- 
ing parameters for spectral modeling and how the spectral fits 
are performed is of some practical importance. We note that 
the most obvious strategy of spectral modeling, namely using 
the best fit values of nearby observations, can severely bias the 
correlations found between different spectral parameters as this 
strategy can force larger numbers of data points onto local x^- 
minima. This bias is less of an issue when always using the 
same starting parameters, which is the approach adopted for 
the fits described below, even though this approach requires a 
larger amount of CPU time. 

We also find that for complex models, such as the 
Comptonization fits described in Sect. 13.41 and 13.51 the x^- 
valley exhibits a large number of local minima and there are 
many degeneracies between parameter combinations. Both ef- 
fects result in a significant fraction of fits converging on lo- 
cal ;i'^-minima. Even when employing standard tricks, such as 

' Tables 1-3 are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
_http://cdsweb.u-strasbg . fr/cgi-bin/qcat?J/A+A/_ 
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computing error contours, the global minimum was not always 
reached. These "rogue fits" can be identified through outliers in 
one or more parameters when looking at correlations between 
different spectral parameters. Where we identified such bad fits 
we refit the data using a different set of starting parameters ob- 
tained from "typical fits" using the same model. 

A potential drawback of this approach is that it effectively 
forces these observations onto the correlations. Only about 
20% of our observations, however, had to be treated in this la- 
bor intensive way. For a small fraction of these refitted observa- 
tions, the x^f.^ of the finally adopted spectral model was found 
to be larger by (9(0.01) than the original due to degenera- 
cies between the fitting parameters. A typical example here are 
fits where the Fe Ka line energy pegs at its (unphysical) lower 
bound of 6keV, which is often the case because of the PC A 
calibration uncertainty around the Xe L-edge. In such cases we 
searched for better models with the Fe line energy at ~6.4keV 
and accepted the best fit parameters from these fits, even if their 

was slightly higher The number of observations for which 
such a strategy was required is still small enough, and the in- 
crease of;^'^ was only (9(0.01), such that this strategy does not 
influence the results listed below and we are confident that all 
parameters shown in the following reflect the physical behavior 
ofCyg X-1. 

In general, for all spectral fits presented here we find that 
for the reduced J < 2. Values of clearly larger than 1 
might appear worrisome at first, however, we stress that the 
signal to noise ratio of our observations is very high and that 
the available spectral models are often too simple to be able to 
describe all subtleties of the observations. Furthermore, given 
our choice of systematics, for lower x^^^ PCA data are com- 
pletely dominated by calibration uncertainties. In principle this 
would require data analysis methods that treat systematic errors 
in a much more elaborate way. We stress that typical observa- 
tions of bright sources with other satellites, such as Chandra 
or XMM-Newton, result in good;^^^ values despite having sig- 
nificantly larger (10% or more) deviations between the data 
and the model than our RXTE fits, where we typically find 
that data and model deviate by <1%. These better values 
from Chandra and XMM-Newton are by virtue of their smaller 
effective areas and consequently larger Poisson errors which 
dominates over the calibration uncertainty of these satellites. 
Seen in this light, our RXTE models thus provide an overall 
good description of the spectral shape of Cygnus X-1. 

3. Modeling the 3-200 keV Spectrum of Cygnus X-1 

3.1. Cygnus X-1 since the launch of RXTE 

Fig lDshows the X-ray and radio lightcurves of Cyg X-1 for the 
nine year long interval from early 1996 (launch of RXTE) until 
the end of 2004. Shown are the 2- lOkeV X-ray flux m easured 
with the RXTE All Sky Monitor jLevine et al. ''l996). as well 
as the radio lightcurves at 15 GHz from the Ryle telescope and 
at 2.25 GHz and 8.3 GHz measured with the Green Bank inter- 
ferometer (GBI; these data are only available until early 2000 
when the GBI monitoring campaign was discontinued due to 
lack of funding). To minimize the significant source variability 



on its orbital time scale jBrocksopp et aljri999ah and to em- 
phasize the long-term evolution of the source, the data were 
rebinned to a resolution of 5.6 d. 

The figure clearly shows the well known fact that Cyg X-1 
transits between two states. As will be elaborated upon below, 
for ASM count rates below ~45 cps, the source is in a classi- 
cal hard state. During this time radio emission is observed at 
a level of ~12mJy at 15 GHz. In contrast, during times when 
the ASM count rate is >80cps, the source spectrum is soft, al- 
though a hard power law component is still seen. We will call 
these phases the "soft state", although it is important to note 
that Cyg X-1 only rarely reaches what one would call a classi- 
cal "soft state" in transient BHC. During these soft states, the 
radio emission is strongly reduced. We especially note that dur- 
ing the long soft state of 2001/2002 there were phases where no 
radio emission could be detected (mm-emiss ion was quenched 
during that time as welklTi gelaar et alJ2004l) . 

During transitions between the hard and soft states, radio 
flares are observed where the radio flux can be up to a factor of 
two higher than during the hard state. Cyg X-1 also shows short 
episodes where the ASM count rate increases. We have iden- 
tified these episodes with "failed state transitions" before (pa- 
per i), since they are characterized by the X-ray spectrum soft- 
ening and by clearly changed timing properties (e.g., increased 
X-ray time lags Pottschmidt et al. 2000), similar to state tran- 
sitions, although the source never really settles into a soft state 
like behavior 

Prior to ~1999, Cyg X-1 was in the hard state for most 
of the time, except for rather short "failed state transitions" 
and the short soft state of 1996. From that time onwards, the 
frequency of the failed state transitions increased until, from 
2001 July onwai-ds, Cyg X-1 reached the long 2001/2002 soft 
state. It stayed in this phase for almost a year until 2002 
August/September where it rapidly transited back into the hard 
state, with a spectrum similar to before the soft state. Shortly 
after this, in 2003 June and at the turn of 2003/2004, two other 
short soft state episodes occured. From early 2004 onwards, the 
RXTE ASM count rate was generally increased with respect to 
the 1997-1999 hard states and the source spectrum was gener- 
ally softer as well. 

We can quantify this behavior for the moment by defin- 
ing the hard state to have an orbit averaged ASM count rate 
of <45 cps (see Sect. 13. 3> and the soft state to have an ASM 
count rate >80cps. Using these definitions, for the total time 
interval from 1996-2004, Cyg X-1 spent 75% of the time in 
the hard state, 4% in the soft state and 21% in between. Before 
MJD51300, the source was in the hard state for 90% of the 
time, in the soft state 4% of the time, and in between for ~6% of 
the time. These numbers change dramatically after MJD 5 1 300, 
where the time spent in the hard state decreased to 66%, while 
the soft and intermediate phases account for 34% of all ASM 
data points. While these numbers are not too precise because 
of our arbitrary cuts between the different states in terms of 
the ASM count rate, they clearly point towards a change in 
source behavior in recent years, which might have been fore- 
shadowed by a change in the timing properties of the source 
around MJD 50920 (see Fig. 3 of paper i). Specifically, around 
this time the short time scale rms variability (i.e., the rms as- 
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Fig.l. RXTE-ASM 2-lOkeV, and ra- 
dio lightcurves of Cyg X-1 from 1996 
until the end of 2004. The 15 GHz 
data are from the Ryle telescope, 
the 8.3 GHz and 2.25 GHz measure- 
ments come from the National Radio 
Astronomy Observatory's Green Bank 
Interferometer. The data are rebinned to 
a resolution of 5.6 d to smooth out the 
orbital variability. 



sociated with individual pointed observations) significantly de- 
creased, while the long time scale variabilty (i.e., the preva- 
lence of failed state transitions) increased. 

The physical reason for this change of behavior after ~MJD 
51000 is unclear. In its simplest interpretation, the soft state 
occurs at higher M, and thus one could speculate that the 
mass accretion rate increased around that date. As M is re- 
lated to the mass loss rate of HDE 226868, in principle it is 
possible to measure M from it s correlation w ith the equiva- 
lent width (EW) of the Ha hne (^ Puls et al.h 996). The simplest 
picture, that a higher mass loss rate of HDE 226868 results 
in a higher M for the black hole, however, is contradicted by 
the observations. Contrary to expectations, th e soft state Ha 



EW i s generally lower than in the hard state (j Voloshinae^al 
'l997l l^ocksopp et al. 1999at iGies etaDl2003t iTarasov et al 



2003h . Fo r single stars, suc h a behavior indicates a lower mass 
loss rate (lGiesetalJl200i and therein). A possible solution 
to this enigma is that photoioniza tion by the X-ray sourc e 
changes the wind ionization state (IHatchett & McCravIl 19771) . 
Observationally, this effect could be the ca use of the varia- 
tion of the Ha hne p rofile with orbital phase dGies et alJl2003t 
Tarasov et al.l 120031) . To explain the relation between the X- 
rays and the optical data, both teams of authors suggest that 
a decrease in the mass loss rate of HDE 226868, M„, changes 
the properties of the accretion flow onto the black hole. A de- 
creased wind density results in an increase of the size of the 
Stromgren sphere around the X-ray source, resulting in a lower 
radiative acceleration and consequential lower terminal wind 
velocity, Vw. Since th e mass accretion r ate is M oc M^/v^, 
(iBondi & Hovlell944tlDavidson & OstrikeJl973l) . an increase 



of M through the accretion disk is triggered, leading to the in- 
crease in soft X-ray emission. The sparse optical coverage pro- 
hibits so far direct tests whether this picture is true, although at 
least for one failed state transition. .Tarasov et al. (.2003) find a 
decrease of the EW before flie flare, consistent with the above 
model. If this general idea is correct, it then follows that the 
recent increase in soft state activity of Cyg X-1 coul d be linked 
to epi sodes of lower mass loss from HDE 226868 dGies et alJ 
I2OO3I) . e.g., linked to activity cycles of the donor star 



3.2. Choice of Spectral Models 

A major issue with describing Galactic black hole X-ray spec- 
tra is the choice of the models to be employed; a multitude of 
more or less physically motivated models are available. It is of- 
ten found that several of these models d escrib e the X-ray data 
equally well in a sense (Nowak et al ]l2002h . Unfortunately, 
however, it can be difficult to directly compare spectral pa- 
rameters obtained from these physically motivated models, as 
they o ften make different implicit assum ptions in their physical 
setup ( Nowak et al. '2002";' 'Coppi'2004' and references therein). 
For example, Comptonization models differ in their assump- 
tions on the spectral shape of the seed photons (e.g., Wien spec- 
trum versus black body versus accretion disk spectrum) and 
they also assume different geometries for th e source of the seed 
photons and of the Comptonizing medium jNowak et aI.l20o3: 
Coppi 2004). These assumptions greatly complicate a direct 
comparison even of "simple" parameters such as the covering 
factor of the reflecting medium. 
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Fig. 2. Unfolded spectrum and residu- 
als for the Cyg X-1 observation of 1999 
July 30 (P40090/14), using the broken 
power law model as the baseline con- 
tinuum. 



Here, we use an approach already u sed by us earlier in a 
study of GX 339-4 dNowak et alJl2002l) to describe the spec- 
tral evolution of Cyg X-1 using several different spectral mod- 
els. We start in Sect. l3.3l with the most simple spectral model, 
a broken power law with an exponential cutoff. This model 
gives a very good empirical description of the hard state X- 
ray spectral shape. We then use two different physically moti- 
vated models to find an interpretation of the empirically derived 
spectral shape in terms of Comptonization. In Sect. 13.41 we 
use th e widel y available Co mptoniz ation rn odel of Titarchuk 
(Il994i see also lTitarchukfe Lvubarsk iil 1 995t lTitarchuk & Hua 
Il995l) . compTT. We contrast the compTT results in Sect. |33] 
with the mo re elaborate Comptonization model eqpair ( Coppi. 
Il992ill999l) . which provides a good example for todays class 
of self consistent Comptonization models that also include the 
capability to generate the spectra resulting from non-thermal 
Comptonization. 

3.3. Empirical models: Brol<en Power Law Fits 

To describe the spectrum of Cyg X-1 on purely empirical 
grounds, as originally motivated by our modeling of the broad 
band radio to X-ra y spectrum of GX 339-4 and Cyg X-1, 
dNowak et alJl2005h . we model the joint PCA/HEXTE data 
with an absorbed broken power law with break energy fibreak, 
soft photon index Fi, and hard photon index that is expo- 
nentially cutoff above an energy Ecm with an e-folding energy 
of Efoid- In addition to this continuum, all fits require an Fe 
Ka line at ~6.4keV, which we model as a Gaussian with en- 
ergy /speKt? and width crpeKa- Table 1, available in electronic 
form only (see footnote [!), contains all of these best fit pa- 
rameters and their 90% confidence intervals, including Fi, F2, 
iifoid, £bieak, and the Fe line parameters. The table also includes 
photon and energy fluxes determined from the model for the 
2-5 keV, 5-lOkeV, 10-50keV, 50-100 keV bands and the un- 



absorbed 2-100 keV flux. Furthermore, the table includes the 
time of each observation, its exposure time, and the PCA and 
HEXTE count rates. 

Fig.|2]shows a typical example for a broken power law fit, 
illustrating that for most observations this continuum model 
gives a good description of the data, with x^^^ < 2. A simi- 
lar result for broken power law fits has been found previously, 
e.g., by Gierliriski et al. ( 1999). In general, we find that dur- 
ing the hard state the goodness of the broken power law fits as 
measured in terms of their x^^^ surpasses that of the compTT 
and eqpair fits. During the soft state, the broken power law 
fits still provide a very good description of the hard spectrum 
above >6 keV, but fare less well below ~6 keV where the ac- 
cretion disk becomes important^. For these data, compTT and 
eqpair give a better description of the spectrum. 

As shown in Fig.|3^, for all observations of Cyg X-1 studied 
here, F2 < Fi, that is, the higher energy power law is spec- 

- This failure of the simple broken power law is most evident for ob- 
servations P501 10/53 (2002.05.06:00), P60090/05 (2002.05.06:02), 
P60090/07 (2002.05.31:23) and P60090/15 (2002.09.20:16). These 
data are characterized by Fi > 3, i.e., they are amongst the softest 
in our sample and there is a clear thermal component present. As our 
main interest in this section is on the behavior of the harder spec- 
trum, we ignored the data below 5keV for these two observations. 
Above that energy, the spectrum could again be well described by a 
broken power law. Additionally, there are two peculiar observations, 
P501 10/07 (2000.06.02:05) and P501 10/10 (2000.07.15:00). In these 
observations, Fi ~ r2 (but not equal!) and when both power law in- 
dices and the Fe line parameters are left free, the uncertainties of both 
power law indices are very large, the break energy in these observa- 
tions converges to ~5.5keV and the energy of the Gaussian is found 
close to 5.5 keV with a width in excess of 1 keV, effectively smoothing 
out the effect of the break. Rather than forcing the power law indices 
to some (arbitrary) value, we have decided to omit these two obser- 
vations from the following analysis, although they are contained in 
Table 1. 
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Fig. 3. a) Correlation between the soft power law index, Fi , and the hard power law index, for the broken power law fits. The 
dashed line corresponds to Fi - Filled diamonds designate hard state observations, defined by Fi < 2.1, open diamonds are 
intermediate and soft state observations, b) Correlation between soft energy power law index, Fi, and power law break, Fi - F2, 
for the broken power law fits. Lines parallel to the dotted line corrspond to lines of constant F2. For both correlations. Spearman's 
rank correlation coefficient p = 0.916. 
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Fig. 4. Correlation between the photon index of the soft power 
law index, Fi , and the e-folding energy of the hard continuum, 
£foid, for the broken power law fits. Symbols are the same as in 
Fig.E] 

trally harder than the soft one. For the majority of our fits, 
the break energy of the power laws is found between 15 and 
40keV. There is a weak correlation between iibreak and fifoid 



in the sense that higher break energies also have higher fold- 
ing energies. A similar correlation exists between Fi and Efoia 
(Fig. |4}, where we find that a softer continuum implies more 
curvature at high energies than harder spectra. In Fig.l^and all 
following figures showing spectral parameters, we show data 
points from the hard state with the filled symbolds and data 
from observations from the intermediate and soft state using 
open symbols. To define the hard state, we use Fi < 2 .1 (a 
simplified version of the classification of Remilla r'dl2005l) . cor- 
responding to ASM count rates below 45 counts s"' . 

The most interesting property of the broken power law fits 
is the hardening of the broken power law, which can be physi- 
cally interpreted as the hardening of the underlying continuum 
caused by Compton reflection off cold or mildly ionized ma- 
terial. As has been pointed out by Zdziarski et al. ( 1999), for 
Seyfert galaxies and galactic black hole binaries, there is a cor- 
relation between the power law index and the covering factor 
of the reflecting medium, Q./2n, in the sense that softer power 
laws tend to show stronger reflection. Our broken power law 
fits show a similar correlation: There is a strong linear relation- 
ship between Fi and F2 (Fig. |3^). Correlating the soft power 
law index, Fi with the hardening of the broken power law at the 
break energy, AF = Fi - F2, results in a simi larly tight cor rela- 
tion (Fig.l^J)). As discussed, e.g., by Zdziars ki et alJ ( fl999t and 
references therein), for standard power law plus reflection fits, 
F and Q/27r are not strictly independent fit parameters: During 

minimization, softer power law continua can be compen- 
sated to some extent by increasing Q.I2n. This behavior can 
lead to artificial correlations between these two parameters. 
While the range of Q and F used to derive the observed Q.- 
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r correlation was much larger than the system atic uncertainty 
caused by the non-independent fit parameters dZdziarski et alJ 
Il999 ). for sources where Q and F do not vary as much, the ar- 
tificial correlation is still of some concern. In contrast, in the 
broken power-law fits Fi and F2 are virtually independent and 
therefore artificial correlations are not an issue for our analysis. 

Lastly, we comment on the behavior of the Fe Ka line 
in the broken power law fits. In general, we find the line 
at 6.4 keV, i.e., consistent with neutral iron, although spectra 
showing a larger AF tend to show slightly higher line ener- 
gies. In no case does the line energy exceed 6.7 keV. This trend 
could indicate a slightly higher degree of disk ionization for 
the softer spectra during these episodes. The distribution of the 
line width, crpe Ka, from all observations can be well described 
by a Gaussian distribution with a mean of 0.6 keV, consistent, 
e.g., with earlier BeppoSAX and RXTE results ( .Frontera et al. 
l2QQl[ ICiilfanovet aljll999l) . We performed Monte Carlo sim- 
ulations, confirming our experience that with its 18% FWHM 
energy resolution at 5.9 keV the PCA can resolve lines down to 
cr ~ 0.3 keV. The fits indicate, therefore, that the Fe Ka line is 
slightly broadened, as has been found also with higher resolu- 
tion instruments (e.g.. Miller et al. 2002) . Note, however, that 
the ~5 keV calibration feature discussed in AppendixElas well 
as the presence of a soft excess can bias the line parameters to- 
wards a lower energy and larger cr. 

3.4. Simple Comptonization Models: compTT 

We now turn to modeling the data with more physically 
motivated spec tral models. We first use the Comptonization 
model compTT ('Titarchuk"l994':'Titar chuk & Lvubarskiill995t 
[Titarchuk & Hua 1995), as historically it has met with good 
success in describing black hole spectra. Model parameters dis- 
cussed here, therefore, can be used in comparisons with earlier 
results and other sources (see, e.g., J^ottschmidt et al. 2003a). 
Furthermore, as we shall see below, the behavior of compTT 
parameters is very similar to that of the more sophisticated 
Comptonization models, while being substantially faster to fit. 

We took compTT as the baseline continuum, assuming 
a disk geometry. The parameters of the continuum are the 
electron optical depth, Tg, and the electron temperature kTf^. 
An accretion disk continuum with a r"^^'* temperature pro- 
file and inner disk temperature kT\„ (XSPEC model diskbb; 
iMitsuda et all Il984t iMakishima' et al.. .1986.) was added to 
model the soft emission. We set the temperature of the seed 
photons for Comptonization equal to kTi„, indicating that the 
disk is the source of the seed photons. Note that this ap- 
proach is only approximately justified, since the seed pho- 
ton distribution assumed by compTT is a Wien distribution 
and not a proper disk spectrum. We take Compton reflec- 
tion of the Comptoniz ation continuum into account using the 
Green's functions of Mag dziarz & Zdziarskil d 19951 XSPEC 
model reflect), assuming solar abundances, a neutral reflec- 
tor, an inclination of 40°, and taking into account absorption in 
the interstellar medium and the stellar wind. The behavior of 
the latter is similar to that seen in the eqpair fits, so we will 
defer a discussion of the behavior of the Hydrogen equivalent 



column, A^H, to Sect. l3.5l All best fit parameters, including kT^n, 
T, kTg, Q/27r, the Fe line parameters, and the model normaliza- 
tions, as well as their 90% confidence intervals are available in 
the electronic Table 2 (see footnote[0. The table also includes 
fluxes from the model in the same standard bands as Table 1, 
as well as the total bolometric flux inferred from the model. 
Finally, the table also includes the total bolometric accretion 
disk flux. 

The Fe Ka line was again modeled by a Gaussian, but with 
the Fe Ka line energy fixed at 6.4 keV (fits with the line energy 
free were virtually identical). There is one clear difference be- 
tween the Fe Ka parameters in the broken power law fits and 
the ones found with compTT: Fe Ka line widths are greater by a 
factor of ~1.55 than in the broken power law fits. Inspection of 
the compTT residuals shows the PCA to have a wavy structure 
at low energies around the Fe line, indicating that our choice 
of soft spectral components does not adequately describe the 
details of the shape of the soft continuum. In addition, we see 
the Ka line width increasing almost linearly with the accretion 
disk flux, until it saturates at Fa^k ~ 10"'^ erg cm"' s"' with a 
line width of ~1 keV. Such a behavior is similar to that of a sim- 
ilarly broad fine feature (at ~5.7keV) seen in XTE J1908-H094, 
in the ear ly and late phases of the outburst of this soft X-ray 
transient tao& ii^ et al.ll2004l) . The broader Fe Ka line width 
can thus be seen as an attempt of the ;^'^-minimization proce- 
dure to smooth out the continuum to provide a more power-law 
like soft continuum. This conclusion is verified by the gener- 
ally worse of the compTT fits compared to the broken power 
law ones, although in general the model is still providing a good 
description of the data: For hard state and intermediate state ob- 
servations, < 1.3, for soft state observations, is clearly 
higher, although (except for four cases) still better than 2. Note 
that only in 12 out of the 202 observations is the x^g^ of the 
compTT fits better than that of the broken power law fits. 

Replacing the disk compon ent with a relatiy istic disk spec- 
trum (XSPEC model diskpn. iGTerliriski et al1ir999.) resulted 
in significantly worse fits. Adding a second soft component 
(modeled either by a second disk spectrurn or by optically 
thick Comptonization iGierliiiski et alJ 1 19991) to our baseline 
diskbb model did not result in improving the residuals sig- 
nificantly either. We therefore decided to stick with our base- 
line model, especially since the Comptonization parameters are 
mainly driven by the harder spectrum. 

Fig- E shows the correlation between the optical depth, r, 
and the temperature kTf^ of the Comptonizing plasma. Defining 
the Compton y parameter in the usual fashion. 



4-kT^ 



(1) 



we find that for a large fraction of the data (~7Q% of the 
observations), y ^ 0.5 (Fig. [S] dashed line). That is, dur- 
ing these times the Compton-y of the Comptonizing plasma 
is nearly constant to within our error bars. This behavior is 
a well known property of plasmas held in thermal equiUb- 
rium by balancing the a vailable heating by Compton cooling 
(e.g.. IPove et al.l [ l9 97a b). Note that we set the lower limit 
of T to 0.01 in our fits. For r < 0.1, if y were to remain 
constant at such low optical depths, the plasma would have 
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Fig. 5. Correlation between the electron temperature of the 
Comptonizing plasma, kTf^ and its optical depth, r. The dashed 
line shows the values for which the Compton y parameter 
y = 0.5, the lower limit for r set in the fits was t - 0.01. 
Symbols used are identical to Fig.|3] 

relativistic tempera tures, violating the app roximations used in 
the compTT model jTitarchuk & Hualll995 ). Processes such as 
photon-photon pair production would become important. For 
these conditions, more realistic Comptonization models such 
as eqpair have to be used (Sect. l33l . We still find that y is 
well determined for these observations, with observations with 
kTe > lOOkeV leaving the track of constant y inferred above 
as T pegs at its lower limit. Virtually all of these observations 
were made during soft state episodes and are therefore charac- 
terized by a strong soft excess and a softer hard spectrum than 
the average hard state spectrum: The majority of spectra with 
a hard power law index of r2 > 1.8 in the broken power law 
fits come from observations for which the compTT fits yield 
> lOOkeV. This result is consisten t with earlier stu dies of 
the 1996 soft state of Cygnus X-1 re.g.. lCui etalJl997l) . 

Since the parameter responsible for the shape of the 
Comptonization spectrum is the Compton-y, we now study the 
relationship between the accretion disk and the reflection pa- 
rameters and y (Figs. |6l and 0. In general, the correlations 
found in this way are much tighter than correlations with re- 
spect to kTf. and t, as these two parameters are strongly cou- 
pled. 

We compute the bolometric unabsorbed flux of the accre- 
tion disk, Fdisk, by integrating the diskbb component of the 
best fit model from 1 e V to 50 ke V. The second source of soft 
photons, the Wien seed photon spectrum, is not taken into ac- 
count. Comparison with the soft fluxes derived with eqpair, 
where the seed photons can be taken into account (see Sect. 13. 51 
below), shows that not taking the Wien spectrum into account 
results in the bolometric disk fluxes from compTT to be 73% of 



Fig. 6. Relationship between Compton-y and the unabsorbed 
accretion disk flux for the compTT fits. Symbols used are the 
same as in Fig.|3] 
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Fig. 7. Correlation between y and the covering factor for re- 
flection, Q/27r for the compTT fits, showing that Q. strongly in- 
creases for small y. Symbols used are the same as in Fig.|3l 

the eqpair fluxes. After correcting the compTT fluxes for this 
offset, the standard deviation of the the eqpair and compTT 
ratio is found to 20%, which we take as the systematic uncer- 
tainty of the individual fluxes caused by extrapolating the best 
fit model outside of the energy range used for fitting. This 20% 
systematic uncertainty of individual flux points is much smafler 
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than the overall variation of the disk fluxes found in the data. 
Note that in the following, when quoting Fdisk from compTT fits 
we do not apply any correction factor, such as might be caused 
by ignoring the Wien seed spectrum. 

Fig Elshows how Fdisk depends ony. There are clearly two 
regimes of different spectral behavior recognizable: During the 
hard state (filled symbols), Fdisk is seen to vary by almost an or- 
der of magnitude, while y remains constant. At the same time, 
Q./2n varies only slightly (Fig.Q. This behavior of parameters 
holds for disk fluxes below ~10"** ergcm"^ s"', corresponding 
to a threshold disk lumino sity of 4.8 x lO^^ ergs'^ (assum- 
ing a distance of 2 kpc, see IZi61kowskill2005l and references 
therein). Above this threshold the intermediate or the soft state 
is reached. Here, y is seen to vary strongly, the disk luminosity 
is far less variable than before, and f2/2;r varies strongly. 

3.5. Comptonization: eqpair 

We turn to the most sophisticat ed of the spectral m odels dis- 
cussed here, the eqpair model llCoDDill992ill999l) . It is a hy- 
brid thermal/nonthermal Comptonization code, and it includes 
electron-positron pair production. Given the overall success 
of thermal Comptonization to describe the hard X-ray spec- 
trum described in Sect. 13.41 in the following we assume a fully 
thermal electron distribution. Eqpair computes the temper- 
ature of the Comptonizing medium self-consistently, by bal- 
ancing external heating with Compton cooling of a defined 
seed-p hoton distribution, wh ich in our case is the diskpn 
model jGierliiiski et aljfl999t) . The amount of heating of the 
Comptonizing medium is specified in eqpair through the ra- 
tio of the compactnesses of the Comptonizing medium and the 
seed photon distribution, where the dimensionless com- 

pactness parameter of a medium of characteristic size r is de- 
fined by 

Lcr-T 

e = — - (2) 

where L is the source luminosity, ctj the Thomson cross sec- 
tion, me the electron mass, and c the speed of light. To force 
the seed photon spectrum to be dominated by disk radiation, 
we fix 4 = 10. A fraction of the hard X-rays is scattered back 
onto the disk where it is reflected. The amount of reflection 
is again quantified by the covering factor for reflection, Q/ln, 
contrary to the compTT fits, relativistic smearing is taken into 
account for the reflection component. 

Similar to the compTT fits, we model the soft excess by 
adding a diskbb to the continuum. In contrast, detectors with 
a lower energy threshold than the PCA often have shown the 
soft X-ray spectrum in eqpair modeling to be more complex 
than what we use here. The soft excess has been described 
by optically thick Comptonization of 1 00-300 eV disk pho- 
tons ( Baluciriska-Church et al. 1995), as predict ed for the inner 
regi ons of Shakura & Sunvaev (1973) disks ( Gier liiiski et alJ 
Il997; .Frontera et aL.200k .Pi Salvo et al...20011 Our spectra 
have an ~3keV lower threshold, thus we cannot describe in 
detail this component. Similar to the compTT fits, we find 
that adding a diskbb spectrum to the eqpair continuum re- 
sults in a satisfactory description of the data. We set the inner 



disk temperature to the maximum temperature of the eqpair 
seed photon distribution. Modeling the soft excess with the 
diskbb works better than with the diskpn, as the shape of the 
diskbb model better approximates the shape of a (saturated) 
Comptonization spectrum. 

Absorption in the interstellar medium and in the stel- 

lar wind is taken into accou nt using the model of 

iBatuciriska-Church & McCammonI lll992h . letting the hydro- 
gen equivalent column vary freely. We find a dependenc e of 
A^H with orbital phase (using the ephemeris of Brocksopp et alJ 
rr999bl) . with A^H being higher by a factor ~2 during superior 
conjunction, pointing towards significant absorption in the stel- 
lar wind of HDE 226868 (see Wen et al. 1999 for a similar con- 
clusion b ased on RXTE-ASM data). We also confirm earlier 
results bv'Balu cinska-Church et al] ( 1200( 1 of a secondary max- 
imum of A^H around orbital phase 0.6, which earlier has been 
identified with absorption in the accretion stream. No other fit 
parameters show a significant dependence on orbital phase. 

Using this approach, our fits gave x^^^ values similar to the 
models discussed in Sects. I?31 and l3.4l with only a very weak 
dependence of x^^^ on the source state (softer spectra tend to 
have slightly worse X^ed)- During the hard state, x^^^ ^'^^ 
eqpair fits is generally worse than that of the broken power 
law fits, while it is better during the soft state (we attribute the 
latter to us not including a thermal component in the broken 
power law fits). All best-fit parameters, including 4/4, ^/2n, 
and the Fe line parameters, and X-ray fluxes for the same bands 
as for the broken powerlaw and compTT fits, as well as unab- 
sorbed fluxes for the accretion disk and for the broadband spec- 
trum are contained in Table 3, which is available in electronic 
form only (see footnote^. 

Fig. |8l shows the evolution of energetics of the soft and 
the hard spectral components. Similar to the compTT fits, we 
characterize the soft flux by the unabsorbed bolometric ac- 
cretion disk flux, Fdisk- This flux is computed as the sum of 
the leV-50keV fluxes of the diskpn^. and diskbb spectral 
components of the best fit models. As discussed in Sect. 13.41 
the systematic uncertainty of these fluxes is estimated to be 
20%, much smaller than the variations observed. Not unex- 
pectedly, the disk flux is correlated with the binary orbit av- 
eraged RXTE-ASM count rate. A more careful analysis of 
the lightcurve reveals, however, that there the ASM does not 
always fully track Fusk- For example, Fdisk increases signif- 
icantly during the second and third of the three dominant soft 
flares around MJD 5 1500 and remains weak during the first. We 
will revisit these differences, which might be due to hysteresis 
eff'ects similar to BHC transient outbursts, below. Overall, how- 
ever, the correlation between Fdisk and the ASM count rate is 
good, with the Spearman rank correlation coefficient indicating 
a highly significant correlation (p = 0.85 1). A linear regression 
gives 



^''"'^ = 0.0592 ^^'^ , - 0.809 
counts s ' 



10 erg cm ^ s '~ ~ ' ^ ^ 

where Fasm is the rebinned ASM count rate. Due to the hys- 
teresis effects mentioned above, we estimate the uncertainty of 

^ The eqpair model has a toggle that allows one to calculate the 
flux of solely the seed photon component 
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Fig. 8. Spectral evolution of Cyg X-1 
from 1999 until the end of 2004. Shown 
from top to bottom are the radio flux 
at 15 GHz as measured with the Ryle 
telescope, the compactness ratio, 
and the bolometric flux of the accre- 
tion disk as obtained from the eqpair 
models, and the 2-lOkeV RXTE- ASM 
light curve, rebinned to a resolution of 
5.6 d to smooth out the orbital variabil- 
ity. Error bars are not shown for clarity. 



Fdisk computed from Eq. |3]to ~5 10"''ergcm"' s"'. We note 
that evolution of the soft excess obtained from the compTT fits 
is very similar to that from the eqpair fits, although in gen- 
eral compTT gives a disk flux that is ~1Q% of that found with 
eqpair. 

The hard spectral component is characterized by the 
compactness ratio, ^h/^s and the optical depth t of the 
Comptonizing medium. Large ^h/^s correspond to the hard 
state. As shown in Fig.|9] t increases with 4/4, indicating the 
increased importance of the Comptonizing medium for higher 
{h/^s- Consequently, low ^h/4 correspond to the soft state, such 
that there is an anticorrelation between ^h/4 and the accretion 
disk luminosity (Fig.llO>. This anticorrelation can be expressed 
as a power law dependence of the bolometric accretion disk lu- 
minosity, Ldisk, and 4/4, 



idisk (4/4) 



-0.19 



(4) 



This relationship holds over the ~1 order of magnitude in vari- 
ation of the disk luminosity. Our fits thus indicate very smooth 
relationships between the parameters of the Comptonizing 
medium and the accretion disk, with no jumps being present 
that would indicate that the source behavior suddenly changes 
when the source switches its state. This behavior is also clearly 
seen in Fig. ^2 where we show how the eqpair parameters re- 
late to the continuum shape as described by the broken power 
law models. For large values of 4/4, the photon index sat- 
urates at ~1.6 and the folding energy approaches its canonical 
value of >150keV. Not unexpectedly, for lower 4/4 as the soft 
state is approached, r2 softens and Em^ decreases, resulting in 
a more curved hard X-ray spectrum. 




Fig. 9. Relationship between the compactness ratio of the 
Comptonizing medium, 4/4, and the optical depth t of the 
Comptonizing medium for the eqpair models. Symbols are 
the same as in Fig.|3] 



We now briefly turn to the evolution of the bo lometric flux 
of Cyg X-1. Based on RXTE and CGRO data. IZhang et alJ 
11 19971) claimed that the difference in bolometric luminosity of 
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Fig. 11. Relationship between and the shape of the X-ray spectrum above lOkeV as found from the broken power law fits, 
a) Correlation between ^h/^s and the photon index, b) Correlation between and the folding energy, Zifoid- 




2 4 6 8 



Fig. 10. The bolometric luminosity of the accretion disk, Ldisk, 
computed assuming a distance of d = 2 kpc, shows a power 
law dependence on the compactness ratio of the Comptonizing 
medium, Symbols used are identical to Fig.|3 

Cyg X-1 between the 1996 hard and soft states was < 50%- 
70%, a statement later revisited by IZdziarski et alJ ll2002) . 
Based on a small number of observations, these authors found 
that the bolometric flux of the source was higher by a factor 
of 3-4 in the soft state. The availability of ou r broad band 
fits allows us to confirm IZdziarski et alJ ll2002h 's earlier re- 



sults: During the RXTE campaign, the mean bolometric unab- 
sorbed flux of Cyg X-1 was 4.3 lO^^ergcm"^ s"', with a mini- 
mum flux of 2.4 lO^^ergcm"^ s"' measured on 2003 January 
10 (observation P60090/23.14off), and a maximum flux of 
9.1 lO^^ergcm"^ s"' measured on 2002 March 25 (observa- 
tion P60090/02.14ofF). Assuming a distance of 2 kpc, these val- 
ues correspond to a mean source luminosity of 2.1 10^^ erg s"', 
varying between 1.2 lO-'^ergs"' and 4.4 10^^ erg s"'. These lu- 
minosity values are extremes; in general we find the luminosity 
of Cyg X-1 during the soft states to be approximately twice that 
of the hard state phases. Cyg X-1 thus has a long-term average 
luminosity of (0.01-0.02) L^m , assuming a black hole mass be- 
tween 10 and 15 solar masses jZiolkow ski 2005:' Herrero et alJ 
.1995> . Fig.^lshows unfolded spectra for two of the more ex- 
treme observations of Cyg X-1, illustrating the strong spectral 
variability of the source. 

Finally, we discuss the variability of the reprocessing fea- 
tures, i.e., the covering factor of the reflecting medium, Q/27r, 
and the Fe Ka line. Fig. 1131 shows that Q./2n ~ 0.15 for 
{h/£s ^ 4, indicating moderate amounts of reflection dur- 
ing the hard state. There is a moderate trend for 0/27r to de- 
crease with increasing (and increasing t), although this 
trend is somewhat weaker and less distinct than previously 
claimed anti-c orrelations between reflection fraction and s pec- 
tral hardness dZdziarski et al.l '200.3': Zd ziarski et 111119991 and 
references therein). For our observations, the eqpair and the 
compTT fits show a similar relative behavior; however, Q./2n 
derived from eqpair is larger by ~0.1. This slightly larger fit- 
ted covering factor is likely due to the usage of a relativistic 
smeared reflection continuum in the eqpair fits; reflection in 
the compTT model is artificially reduced to minimize the sharp 
Fe edge in the reflection model used with those latter fits. There 
may also be influences from the different shapes of the expo- 
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Fig. 12. Unfolded spectra and residu- 
als for the eqpair fits of the Cyg X- 
1 observations of 2003 January 10 
(P60090/23.14off; hard spectrum) and 
2003 July 29 (P60090/35.14off; soft 
spectrum), illustrating the typical spec- 
tral variability of the source. 



♦ 
♦ 



_i I i_ 



_i I i_ 



the spectral parameters considered here. This lack of correla- 
tion is puzzling. Since the Fe Ka line is thought to be the result 
of fluorescence in the accretion disk, we would expect the Fe 
line to track Q./2n, i.e., the variation of the Compton reflection 
component. Similar to our fits with the broken power law or the 
compTT model, there is a larger fraction of fits in which the Fe 
Ka line has an energy that is significantly below 6.4 keV. The 
residuals for these observations often show the Xe L edge cal- 
ibration feature, indicating that calibration effects might lead 
at least in part of the observations to spurious Fe Ka line pa- 
rameters, although we cannot exclude other eff'ects, such as a 
more complex shape of the soft excess. As already indicated 
above, however, neithe r the addition of a Compton component 
(Gierliriskiet al l 19991) nor the substitution of the di skbb soft 
component by a di skpn component result in a better descrip- 
tion of the soft excess and the Fe Ka line. 



3.6. Disk-Jet Interaction 



Fig. 13. Correlation between the compactness ratio of the 
Comptonizing medium, and the covering factor of the 

reflector, Q./2n. Symbols are the same as in Fig.|3] 

nential cutoff in these two models. In the eqpair models, as the 
source approaches the soft state and (h/£s decreases, there is no 
clear correlation with Q./2jt. The reflection fraction varies in a 
seemingly independent manner, although in general our high- 
est fitted reflection fractions come from the soft state. Again 
this is reminiscent of the behavior already seen in the compTT 
fits (Fig.|7}. 

The second indicator for reprocessing in cold material, the 
Fe Ka Une, does not show a correlated behavior with any of 



As we have seen in Fig.|8] in agreement with studies of other 
sources, we find evidence for correlated radio and X-ray vari- 
ability. Due to the comparably low luminosity variation of a 
factor of ~4, Cyg X-1 is not an ideal source to study disk-jet 
correlations across widely separated states, as occur in X-ray 
transients that show factors of 100 to 1000 l uminosity changes 
dFender & Bellonil2004lFender et al.l2004 . The proximity of 
Cyg X- 1 to the threshold between the hard and soft states, how- 
ever, allows us to study this crucial transition interval in much 
greater detail than for transients. 

To parameterize our data with the least possible theoret- 
ical bias, we can use the results from the broken power law 
fits. As shown in Sect. 13.31 these fits indicate the presence of 
two spectral components that vary in a correlated manner as 
described by the Fi-Fa correlation. We show in Fig. [l4l that 
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Fig. 14. The power law break, AF, of the broken power law fits 
can be interpreted as a measure of the bolometric accretion disk 
luminosity. 



Fig. 15. For power law breaks <0.5 the radio flux varies seem- 
ingly independent of the X-ray color, expressed here as the ratio 
between the 2-lOkeV and the 10-50keV flux. Symbols used 
are identical to Fig.|3l See text for further explanation. 



there is an almost linear relation between the soft excess lumi- 
nosity found from Comptonization fits and AF - Fi - F2. In 
Comptonization models, therefore, one can interpret the power 
law break in part as a measure of the bolometric accretion disk 
luminosity, although part of the break can also be caused by the 
Q-F-relationship. 

Alternatively, however, in models explaining the X-ray 
spectra of BHC in terms of emission from along a jet, the 
power law break describes the transition from synchrotron ra- 
diation (as well as accretion disk emission) to a synchrotron 
self-Compto n spectrum that is modified by Compto n re- 
flection (.Markoff et aljbool EoOS: Markoff & Nowaki2004t 
iNowak et alJl2005h . Figure 1151 provides an important bench- 
mark for such models. For AF < 0.5, i.e., most hard state obser- 
vations, the radio flux is seemingly independent of AF, while 
above that threshold strong variations of the radio flux are ob- 
served. 

Using the Comptonization interpretation of the X-ray 
spectrum, we can study how the radio emission and the 
Comptonizing medium interact. For such a study the eqpair 
models are especially interesting, since 4/4 is a measure of 
the energetics of the Comptonizing medium. We find that the 
radio flux is highest at intermediate 4/4 (4/4 ~ 3; Fig.ll6>. 
This behavior is a consequence of the "failed state transitions", 
where Cyg X-1 shows characteristic flaring behavior in the ra- 



dio, similar to what has also been seen in GX 339-4 and other 


iet sources JCorbel et al...200Cl 


2003; Hannikainen et al...l998; 


iFender et alJ 


12004 


see also Figs. [T] and |8}. Note fliat 4/4 = 3 



corresponds to a soft photon index of Fi ~ 2.1, consiste nt with 
the up per limit for the hard state Fi as defined by .Remillar j 
(l2005h . Note also that there seems to be a much tighter rela- 
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Fig. 16. Correlation between the 15 GHz radio flux and the 
compactness ratio of the Comptonizing medium, 4/4- The in- 
crease in radio flux for intermediate values of 4/4 is due to 
the radio flares usually associated with "failed state transitions" 
and the intermediate state. The radio flux shown is the mean ra- 
dio flux averaged over the 5.6 d orbital timescale surrounding 
each RXTE measurement. Symbols used are the same as in 
Fig.E] 
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to be sligh t ly diff eren t from that of X-ra y transients, where 
iGallo et alJ ( l2003h and lCorbel et"al] (l2003l) find that the radio 
flux, Fr, and X-ray flux. Fx are related by oc F^^ . Such 
a scaling is expected in various models in which the X-ray 
and radio producing r egions are connected, either by the X- 
rays originating in a jet fralcke & BiermannllQQfitFcorbel etall 
E003; Markoffetal. 2003, 2005) or by the X-rays originat- 
ing in an advection dominat ed type of accretion flo w that is 
coupled to the radio outflow ( Heinz & Sunyaevl2003l) . Fig. 1171 
shows that for Cyg X- 1 there are indications that the radio-X- 
ray relationship is steeper than this canonical value. Taking all 
observations with available radio fluxes into account, we find 
oc f'^^ (Fig. El solid line). Limiting ourselves to the hard 
state observations only, where the radio fluxes are best deter- 
mined, we find oc F^"^ (Fig. El dashed line), which, given 
the uncertainty of the flux determination, can be considered 
as consistent with the Fr oc f relationship. Similarly steep 
radio-X-ray relationships have been seen in the black hole can- 
didate XTE J 1908 +094 (Jonker et al. 2004) and in the neutron 
star system 4U 1728-34 dMigliari et al.l2003l) . 



Fig. 17. Correlation between the 10-50keV hard X-ray flux 
and the 15 GHz radio flux. Th e dotted line corresponds to the 
canonical Fr oc F^ relation dCorbel eta l. 2003; ' Gallo et al.l 
12003 ). the dashed line is Fr oc F^^ determined from the hard 
state data, the solid line is Fr oc F^^ found from all Cyg X- 1 
measurements. 



tionship between the radio flux and for the soft state com- 
pared to the scatter seen for the hard state. 

In general, BHC for which radio fluxes have been measured 
exhi bit a strong correlation between the r adio and the X-ray 
flux dOaUo et alJ2003UMerloni et alJ2003h . For Cyg X-1, this 
correlation become s less we ll pronounced as the soft X-ray flux 
increases (paper iii: lGallo et al. 2003) . We likewise do not find 
a significant relationship between the radio flux and either the 
bolometric flux o r the disk emission. As already described in 
paper iii and by Nowak et al. however, the radio flux 

correlates well with the 10-50keV hard X-ray flux (p = 0.625). 
For 15 GHz fluxes above 0.1 mJy, the lowest flux at which the 
source was detected during the campaign, we find that the 10- 
50 keV X-ray flux, Fio-so, can be found from 



, F 10-50 

logio 

erg cm-^ s ' 



= -8.27 -H 0.320 logi, 



^ ISGHz 

mJy 



(5) 



with an uncertainty of less than 15%. Here, the radio flux is 
again the mean flux determined over one orbital period of the 
systern to average out the orbital variability ( Brocksopp et alJ 
119993*). 

This result is in strong support of the idea that the hard 
X-ray flux is tightly coupled to the radio, suggesting a possible 
physical connection between the regions in which this radiation 
is produced. Note that for many X-ray transients the radio-X- 
ray relationship is found from the 2-lOkeV ]?XTE-ASM data, 
which is possible only because of the less complex soft ex- 
cess in these Roche lobe accreting systems. Overall, however, 
the radio behavior of Cyg X-1 as expressed by Eq.|5l seems 



4. Summary 

In this paper we have analyzed ~200 observations of 
Cygnus X-1 in all of its states. The spectral parameters are 
in qualitative agreement with those found for Comptonization 
models of earlier observations (e.g., Pottschmidt et al. 2003a; 



Pi Salvo et aL^OOl; Frontera et al. 200 1; Gilfanov et al. 1999; 
Gierlihski et aljl999HDove et al.il998h . The availabihty of flie 



large number of observations allows us to study the variation of 
the spectral parameters in this black hole candidate in unprece- 
dented detail. This has revealed a large number of interesting 
correlations that must be explained by models for the emission 
process and emission geometry of Cyg X-1. In summary, the 
main results of the spectral analysis are as follows: 

1. Cyg X-1 changed its behavior in early 2000. While it was 
predominantly a hard state source before that, since 2000 
the source spent ~34% of the time in the intermediate and 
the soft state (see also paper i). The physical reason for this 
change in behavior could be related to changes in the mass 
loss rate from HDE 226868. 

2. The spectral continuum of Cyg X-1 is very simple and can 
be well described at the level of the RXTE resolution by a 
broken power law with an exponential cutoff. 

3. More physical models such as the two Comptonization 
models considered here try to reproduce this simplicity 
and, for the hard state, result in values comparable (but 
slightly worse) than the broken power law fits. 

4. At the;t'^ level, simple and advanced Comptonization mod- 
els cannot be distinguished and the qualitative variation of 
their parameters is similar From our modeling, a purely 
thermal electron distribution seems to be sufficient to de- 
scribe the spectra in the hard and soft state observations, at 
least in the ~3-120keV energy range considered here. 

5. There is a continuum of spectra between the hard state and 
the soft state, with no strong discontinuities found in the 
behavior of the spectral parameters. This continuous vari- 
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ation is in strong contrast with the changes in timing be- 
havior, where the states are clearly distinguishable. For ex- 
ample, the intermediate state and failed state transitions 
have increased X-ray time lags ( Pottschmidt et al. 2000^ 
and changed power spectra (paper i). 

6. This continuous distribution of parameters from the soft to 
the hard state leads to robust correlations between spec- 
tral parameters, reflecting trends in the continuum of spec- 
tral shapes of Cyg X-1. For the broken power law fits, the 
most important correlation is that between the soft and the 
hard power law index, for the Comptonization fits there is a 
power law relationship between the soft excess luminosity 
and the compactness ratio, L^sk (^h/^s) ° 

7. For the compTT models, the Comptonizing plasma shows 
a steady behavior during the hard state and strongly varies 
during the soft state (with a similar behavior in the reflec- 
tion covering factor, Q/27r). During the hard state, the soft 
excess properties vary seemingly independently from the 
Comptonizing plasma. 

8. Although both Comptonization models show a slight anti- 
correlation between spectral hardening and reflection frac- 
tion, Q/ 27r, it is not as pronou nced as has been previously 
claimed ( Zdziarski et aLbOoil and references therein). It is 
likely that part of this correlation is in fact related to corre- 
lations between broad-band soft and hard continuum com- 
ponents. 

9. The Fe Ka line is found to be only moderately broad in 
the broken power law fits, consistent with earlier results 
llFronteraet al.ll200ll: ICiilfanov et aLlll999l) . An increased 
broadening seen in the Comptonization models is likely not 
a physical effect but rather is caused by an attempt of the fits 
to improve the modeling of the soft excess component. 

10. For the hard state, Cyg X-1 follows the general radio-X- 
ray relation, oc f'^'', although there are indications that 
the relationship is slightly steeper when also taking the soft 
state data into account. 

A further interpretation of these results in terms of physical 
models requires the consideration of the simultaneous variation 
of the timing parameters, which we will present in a further 
paper in this series. 
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Appendix A: Revisiting the RXTE spectral 
calibration 

In this Appendix we summarize our results for the spectral cal- 
ibration ofdieR2£TE based on HEASOFT 5.3.1. As shown, 
e.g., by IWilms et alJ {1999'), earlier HEASOFT versions ex- 
hibited a significant difference in spectral slope of ~0.05 be- 
tween the PCA and HEXTE, with PCA spectra being consis- 
tently softer than spectra derived from HEXTE and other mis- 
sions. For observations of black holes, such a difference is es- 
pecially worrisome since the apparent hardening implied by the 
older calibration can mimic a Compton reflection component. 
In HEASOFT 5.3.1, this discrepancy is now much smaller and 
fits to RXTE observations of the Crab nebula and pulsar give 
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Fig.A.l. Ratio between data and the two power-law spectral 
model for the joint PCA-HEXTE fit of the RXTE observation 
of the Crab nebula and pulsar 40805-01-05-01 described in the 
text. For clarity, only the PCA residuals are shown. The error 
bars reflect the 0.5% systematics applied in our analysis. 

the canonical spectral shape (see ljahoda et alJl2005l for an ex- 
tensive discussion of the PCA calibration). For example, mod- 
eling the joint PCA and HEXTE data from RXTE's Crab obser- 
vation 40805-01-05-01 in 1999 we find that the joint spectrum 
can be described as the sum of two power laws with photon 
indices Fi = 1.45 and = 2.12. The nebula flux at IkeV is 
found to be 10.3 photons cm ~^ s~' keV~', or about 8% higher 
than the canonical value of iToor & SewardI (Il974l) . the hard 
pulsar continuum contributes 0.05 photons cm" -^s-'keV"' at 
this energy (all values given are unabsorbed fluxes, absorption 
in the interstellar medium was taken into account by using the 
cross-sections and abundances of Wilms et al. 2000 and by as- 
suming a H-equivalent column of A^h = 4.2 4 x 10^' cm"^ as 
found with Chandra bv lWeisskopf et al.l2004ii 

The ratio between the best fit Crab model (folded through 
the response matrix) and the data reveals the presence of en- 
ergy dependent systematic uncertainties (Fig. lA.U . These are 
strongest in the 4-6 keV band, around the Xe L-edge, where 
the deviation can amount to up to 1%, and above ~25 keV. We 
chose to take these deviations into account by adding an en- 
ergy independent uncertainty of 0.5% to the data. Applying 
this uncertainty to the Crab data gives a reduced chi-squared 
°^^red ~ 1-38. While a larger systematic uncertainty in the 4- 
6keV band and above ~20keV would be justifiable in light of 
this result, we have decided not to include energy dependent 
systematics since experimenting showed that such systematics 
would affect the Fe Ka line region even more strongly than our 
current strategy (where the Fe Ka line energy is already biased 
towards energies below 6.4 keV). In addition, we only use PCA 
data from channels 6-52 of the standard2f data mode (approx- 
imately 3keV-25keV). For the HEXTE we use channels 17- 
117 (corresponding to ~18keV to ~120keV). The upper en- 
ergy bound of 120keV is set by the shorter observations in the 
sample and by the desire to have an uniform energy range for 
all observations. Note that HEXTE residuals in channels 17-20 
typically show a linear increasing trend, which might be due to 
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calibration effects. We chose to include these channels anyway 
to increase the overlap between the PCA and HEXTE. When 
ignoring these HEXTE channels, thex^ of our fits typically de- 
creases without significant changes of the best fit parameters. 

In addition to the difference in spectral slope, earlier 
HEASOFT versions also showed a very significant overesti- 
mate of the PCA flux with respect to the HEXTE, with the PCA 
fluxes being consistently higher by ~3Q%. This overestimate 
can be taken into account in -fitting by multiplying the spec- 
tral model instrument dependent constant. Traditionally, this 
normalization constant was set to 1 for the PCA and let vary for 
the HEXTE. This traditional approach is unfortunate, since it 
implies that the flux normalization of the HEXTE is uncertain, 
while in reality it was the PCA that showed clear deviations 
in flux with respect to other instruments. For consistency with 
other studies, however, here we continue using the traditional 
approach. For all different spectral models used in this paper, 
we find that under HEASOFT 5.3.1 both instruments now agree 
remarkably well, with the distribution of the constant values 
being consistent with a normal distribution of mean 1.01 and 
standard deviation 0.03. We still chose to include the constant 
in our final fits, as in individual observations a mismatch be- 
tween the PCA and the HEXTE is possible, for example for 
those observations where the PCA deadtime becomes impor- 
tant"^. There are only 7 observations with constants more than 
2cr away from the mean. Except for the significantly deviating 
normalization constant, these observations are not remarkable 
in any other way, and are therefore included in our further anal- 
ysis, although the cause for the deviation is currently unknown. 
We conclude, therefore, that for the overwhelming majority of 
observations there is now agreement in the PCA and HEXTE 
derived fluxes. 

Appendix B: Confidence Intervals and RXTE 
Systematic Errors 

As shown in the previous appendix, the high signal to noise 
observations from RXTE require us to take the systematic er- 
ror of the instrument into account in order to obtain reasonable 
spectral parameters. In general, this is done by adding the sys- 
tematic error in quadrature to the Poisson errors estimated from 
the observation. For instruments such as the PCA it can happen 
that the systematic error starts dominating the overall error. As 
a consequence, the standard "x^in AAf^''"Procedure_for ob- 
taining co nfidence intervals, as ou t lined, e.g., bv lLampton et alJ 
lll976h or lBevingto n & Rob insonI (Il992h . does not apply since 
it is based on the assumption that the variance of the data points 
is due to a Gaussian distribution. Since the addition of system- 
atic errors decreases the x^ sum, confidence intervals estimated 
using the Lampton et al. ( 1976) prescription in the presence of 
systematic errors are significantly larger than what is justified 
from statistical grounds. 

PCA spectra were not deadtime corrected as the spectral shape is 
not influenced by the deadtime correction and as good indicators for 
the PCA deadtime are not available for the time since the failure of the 
Xe layer of PCUO, which also changed the deadtime behavior of that 
PCU. 



Table B.l. Error bars for a representat ive observation of O 
X-1, computed using the prescription of lLampton etalJllT97( 
See text for further explanation. 



Parameter 


value 




90% error bar 








with 


without 








systematics 


kT,, [keV] 


0.855 


(- 


-0.02, +0.03) 


(-0.01, +0.06) 


^diskbb 


1276 


(- 


-175, +203) 


(-24, +40) 


kT, [keV] 


108 


(- 


-31, +45) 


(-20, +27) 


T 


0.26 


(- 


-0.14, +0.13) 


(-0.12, +0.15) 


^comptt [10 "] 


1.86 


(- 


-1.0,+1.1) 


(-0.5, +0.7) 


n/2n [%] 


17.9 


(- 


-1.8,+1.3) 


(-1.1, +1.1) 


CTpcKn [keV] 


1.01 


(- 


-0.05, +0.04) 


(-0.006, +0.012) 


ApeK. [10-2] 


5.34 


(- 


-0.4, +0.4) 


(-0.07, +0.08) 


Chexte 


0.999 


(- 


-0.004, +0.005) 


(-0.003, +0.003) 



To study the influence of the systematic error on the con- 
fidence intervals in our observations, we have used a Monte 
Carlo simulation: We first determined the best fit parameters for 
the comp TT model of observat ion P60090/16 (2002.10.06:15) 
using the lLampton et alJlll976h approach, including systematic 
errors. We then used this best fit model to simulate how this ob- 
servation would look like if the PCA and HEXTE calibration 
were perfectly known. These simulated spectra were then refit 
without applying a systematic error and new confidence inter- 
vals were determined. The resulting 90% confidence intervals 
for both fits are shown in Table lB.l1 (Error bars at the 68% level 
can be obtained from the 90% confidence intervals to a suf- 
ficiently high precision by assuming that the uncertainty can 
be approximated by a normal distribution). Depending on the 
parameter, the error bars are smaller by a factor of 2 to 5, show- 
ing that the confidence interval is completely dominated by the 
systematic error 

Since the systematic error will affect all observations of the 
Cyg X-1 RXTE campaign in a similar way and since the ma- 
jor interest is in determining trends in the evolution of spec- 
tral parameters, using the overestimation of the range of the 
confidence intervals determined from individual observations 
would lead to misleading results. It is better, therefore, to talk 
about an "ensemble average" in determining the errors and to 
take the scatter found in individual relationships between the 
fit parameters of Cyg X-1 in a similar state as an estimate their 
uncertainty. For this reason, we have decided to not show the 
misleading individual error bars in the figures in the main part 
of this paper, although they are available in the online data ac- 
companying this paper (Tables 1-3). 



